arXiv:hep-th/0512108vl 9 Dec 2005 


Application of the Renormalization-group Method to the Reduction 

of Transport Equations 

Teiji Kunihiro 1 and Kyosuke Tsumura 2 

1 Yukawa Institute for Theoretical Physics, Kyoto University, Sakyo-ku, Kyoto 606-8502, Japan 
department of Physics, Kyoto University, Sakyo-ku, Kyoto 606-8502, Japan 

Abstract 

We first give a comprehensive review of the renormalization group method for global 
and asymptotic analysis, putting an emphasis on the relevance to the classical theory of 
envelopes and on the importance of the existence of invariant manifolds of the dynamics under 
consideration. We clarify that an essential point of the method is to convert the problem from 
solving differential equations to obtaining suitable initial (or boundary) conditions: The RG 
equation determines the slow motion of the would-be integral constants in the unperturbative 
solution on the invariant manifold. The RG method is applied to derive the Navier-Stokes 
equation from the Boltzmann equation, as an example of the reduction of dynamics. We 
work out to obtain the transport coefficients in terms of the one-body distribution function. 


1 Introduction 

The concept of the RG was introduced by Stuckelberg and Petermann as well as Gell-Mann and 
Low[I] in relation to an ambiguity in the renomalization procedure of the perturbation series in 
quantum held theory (QFT). However, the essential nature of the RG is exact and hence non- 
perturbative, which was revealed and emphasized by Wilson[2]. Subsequently, as is well known, 
the machinery of the RG has been applied to various problems in QFT and statistical physics 
with a great success. Jlj 

The essence of the RG in quantum held theory (QFT) and statistical physics may be stated 
as follows: Let T(0, gr(A),A) be the effective action (or thermodynamical potential) obtained by 
integration of the held variable with the energy scale down to A from infinity or a very large cutoff 
Ao- Here g( A) is a collection of the coupling constants including the wave-function renormalization 
constant dehned at the energy scale at A. Then the RG equation may be expressed as a simple 
fact that the effective action as a functional of the held variable 0 should be the same, irrespective 
to how much the integration of the held variable is achieved, i.e., 

r(0, 9 (A),A) = r(0,g(A'),A'). ( 1 . 1 ) 

If we take the limit A' —> A, we have 


dT((f),g( A), A) _ Q (1 . ?) 

dA 

which is the Wilson RG equation|2|, or the how equation in the Wegner’s terminology [3j; notice 
that Eq. m is rewritten as 


(K dg _ (W 
dg dA dA 


(1.3) 


If the number of the coupling constants is finite, the theory is called renormalizable. In this case, 
the functional space of the theory does not change in the how given by the variation of A. 
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Owing to the very non-perturbative nature, the RG has at least two merits: (A) Resummation 
of the perturbation series Applying the RG equation of Gell-Mann-Low type [I 3 to pertur¬ 
bative calculations up to first lowest orders, a resummation in the infinite order of diagrams of 
some kind can be achieved. That is, the RG method gives a powerful resummation method [I]. 
(B) Construction of infrared effective actions The RG of Wilson typejEj provides us with a 
systematic method for constructing low-energy effective actions which are asymptotically valid in 
the low-frequency and long-wave length limit. 

An appearance of diverging series is a common phenomenon in every mathematical sciences 
not restricted in QFT, and some convenient “resummation” methods are need and develop ed@. 
Deducing a slow and long-wave length motion is one of the basic problems in almost all the fields 
of physics; for example, statistical physics including the physics of the pattern formation and the 
theories of collective motion in a many body system. The problems may be collectively called the 
reduction problem of dynamics. The RG method^ 0 111 13113113 U2j might be a unified method 
for the reduction of dynamics as well as a powerful resummation method. 





Figure 1: The geometrical image of the reduction of the dynamics. The dynamical variable W (t) 
in the n-dimensional phase space approaches to and after some time is eventually confined in the 
well-defined manifold M as f increases. 

It is noteworthy that one can draw a clear geometrical image for the reduction of dynamical 
systems. Let W(f) be an n-dimensional dynamical system governed by the evolution equation, 

dW 

= F(W, t), (1.4) 

where F is an n-dimensional vector; the dimension n may be finite or infinite. When the dynamics 
is reduced to an m-dimensional system with m being smaller than or equal to n, the vector W(f) 
approaches a well-defined m-dimensional manifold M embedded in the n-dimensional phase space, 
as shown in Fig. 1; then the geometrical object M is called an attractive manifold. If after some 
time the W(f) is confined in the manifold M, M is called an invariant manifold. Furthermore, 
when the dynamics on M is slow, M is also called a slow manifold. Let any point W within the 
manifold M be given by the relation W = R(s) with an m-dimensional parameter s(t). Then the 
reduced dynamics is given by 

d s 

- = G( S ), W = R(s), (1.5) 

where the first equation with the vector field G(s ) defined on M gives the reduced dynamics 
within the manifold M and the second is the representation of M. In the quantum field theory, 
the dynamical variable W corresponds to the set of coupling constants and the renomalizability 
may be interpreted as the existence of an invariant manifold in the space of coupling constants; 
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the increasing time t and the vector held G(s) correspond to the decreasing energy cutoff A and 
the (3 function, respectively. 

In this article, (1) we show that the RG gives a powerful and systematic method for the 
reduction of dynamics and also provides a transparent way for the construction of the attractive 
slow manifold. (2) we apply the method to have the fluid dynamical limit of the Boltzmann 
equation as an example of the reduction of dynamics and the construction of the slow manifold [12]. 
We derive the Navier-Stokes equation explicitly from the Boltzmann equation for the first time; 
the microscopic expressions of the transport coefficients are given. (3) We will put an emphasis 
on the relation of the underlying mathematics of the RG method with the classical theory of 
envelopes in mathematical analysis 0 ei ei mu eh. 


2 The RG method and the classical theory of envelopes 


We here give a brief review of the theory of envelopes. Although the theory can be formulated in 
higher dimensions 0 JH, we consider here envelope curves, for simplicity. 

Let {C T } T be a family of curves parameterized by r in the x-y plane; here C T is represented by 
the equation F(x,y,r) = 0. We suppose that {C T } T has the envelope E , which is represented by 
the equation G(x,y) = 0, as shown in Fig.2. The problem is to obtain G(x,y ) from F(x,y,r). 

Now let E and a curve C TQ have the common tangent line at (x,y) = (xo,yo), he., (xo,2/o) 
is the point of tangency. Then xo and yo are functions of ro; xq = </>(to), yo = V ; ( r o)i an d 
of course G(x 0 ,y 0 ) = 0. Conversely, for each point (x 0 ,y 0 ) on E, there exists a parameter r 0 . 
So we can reduce the problem to get r 0 as a function of (x 0 ,y 0 ); then G(x,y ) is obtained as 
F(x, y, t(x, y)) = G(x,y). Notice that since there is a relation G(x 0 ,y 0 ) = 0 between x 0 and y 0 , 
To is actually a function of Xo or y 0 . To(xo,yo ) can be obtained as follows. 

Since the tangent line of E at (xo, yo) is perpendicular to the normal direction of F(x, y, r) = 0 at 
the same point, one has F x (x o, yo, To)4 > '{ T o) + F y (xo, yo, To)'0 / ('To) = 0. On the other hand, differenti¬ 
ating F(x(r 0 ), y(r 0 ), r 0 ) = 0 with respect to r 0 , one also has F x (x 0 , y 0 , r Q )(f) , (To)+F y (xo, y 0 , r 0 )^'(To)+ 
F tq (xq, y 0 , r o) — 0. Combining the last two equations, we have 


Fro 0&O, yo, A)) 


<9F(x q ,?/o,to) 

8to 


( 2 . 1 ) 


This is the basic equation of the theory of envelopes. Notice that the envelope equation has the 
similar form as the RG equation. 



Figure 2: A family of curves F(x,y,To) = 0 parameterized with To and its envelope defined by 
G(x,y) = 0. 
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One can thus eliminate the parameter r 0 to get a relation between x 0 and y 0 ; G(x,y) = 
F(x,y,r 0 (x,y)) = 0, with the replacement (x 0 ,y 0 ) —> (x,y). G(x,y ) is called the discriminant of 
F(x,y,t). 

When the function F has an additional dependence on a vector C(r), i.e., F = F(x, y, r, C(r)), 
the envelope equation reads 


F TO (x 0 ,yo,T 0 ,C(T 0 )) 


dF(x 0 ,y 0 ,To) dC <9F(x 0 , ?/ 0 , r 0 , C(t 0 )) 
<9t 0 <9t 0 dC 


( 2 . 2 ) 


Comments are in order here: (1) When the family of curves is given by the function y = f(x, r), 
the envelope equation is reduced to df /dr ,o = 0; the envelope is given by y — f(x, To(x)). (2) The 
equation G(x, y) — 0 may give not only the envelope E but also a set of singularities of the curves 

{C T } T . 


3 The RG method; a simplest example 


In this section, using a simplest example we show how the RG method works for obtaining global 
and asymptotic behavior of solutions of differential equations. We shall present the method so 
that the reader will readily see that the notion of envelopes is intrinsically related to the method. 
We shall emphasize that an essential point of the method is tuning the initial condition at an 
arbitrary time to perturbatively along with solving the perturbative equations successively. One 
will see that the reasoning for various steps in the procedure and the underlying picture are quite 
different from the original ones given in jHj. We believe, however, that the present formulation 
emphasizing the role of initial conditions and the relevance to envelopes of perturbative local 
solutions straightens the original argument, and is the most comprehensive one. 

Let us take the following simplest example to show our method: 


d 2 x 

dt 2 


dx 

e— 

dt 


+ x 


0, 


(3.1) 


where e is supposed to be small. The solution to Eq. CB reads x(t) = Rexp(—|t) sin(^/l — ^ t+6 ), 
where A and 6 are constants. 

Now, let us obtain the solution around the initial time t = to in a perturbative way, expanding 
x as x(t,t 0 ) = xo(t,to) + exi(t,to) + e 2 X 2 (t,t 0 ) + where x n (t,t 0 ) (n = 0,1,2...) satisfy 
%o + x o — 0, X n+ \ + x n+ \ = — x n . 

The initial condition may be specified by 


x{t 0 ,t 0 ) = W(t 0 ). (3.2) 

We suppose that the initial value W(t 0 ) is always on an exact solution of Eq. (1.3. If) for any t 0 . We 
also expand the initial value W(t 0 ); VE(t 0 ) = W 0 (to) + eWi(t 0 ) + e 2 W 2 {t 0 ) + ..., and the terms 
Wi(to) will be determined so that the perturbative solutions around different initial times to have 
an envelope. Hence the initial value W(t) thus constructed will give the (approximate but) global 
solution of the equation. 

Let us perform the above program. The lowest solution may be given by x 0 (t,to) — A(t 0 ) sin(t + 
6(to)), where we have made it explicit that the constants A and 6 may depend on the initial time to- 
The initial value W(to) as a function of t 0 is specified as Wo(to) = xo(to, to) = A(t 0 ) sin(t 0 + d(t 0 )). 
We remark that the zero-th order solution a neutrally stable solution; with the perturbation 
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e 7 ^ 0 the constants A and 9 may move slowly. We shall see that the envelope equation gives the 
equations describing the slow motion of A and 9. 

The first order equation now reads x,\ + X\ = —Acos(t + 9), and we choose the solution in 
the following form, x\(t,to) = -4-(t — to) sin(t + 9), which means that the first order initial 
value W\{to) = 0 so that the lowest order value Wo (to) approximates the exact value as closely as 
possible. Similarly, the second order solution may be given by X 2 (t) = j{(t — to) 2 sin (t + 9) — (t — 
t 0 ) cos(t + 0 )}, thus W 2 (to) = 0 again for the present linear equation. 

It should be noted that the secular terms have appeared in the higher order terms, which are 
absent in the exact solution and invalidates the perturbation theory for t far from to- However, with 
the very existence of the secular terms, we could make Wi(tp) (i = 1, 2 ) vanish and W (to) = Wo (to) 
up to the third order. 

Collecting the terms, we have 

A 

x(t,to) = Asin(f + 9) — e— (t — to) sin(t + 9) 

A 

+e 2 — {(t - t 0 ) 2 sin (t + 9) - (t - t 0 ) cos (t + 9)}, (3.3) 

O 


and more importantly W(tp) = Wo (to) = A(to) sin(t 0 + 9 (to)), up to 0(e 3 ). We remark that 
W(to) describing the solution is parameterized by possibly slowly moving variable A(tp) and 
0 (t o ) = t 0 + 9(t 0 ) in a definite way. 

Now we have a family of curves {C to } to given by functions {x(t,t 0 )}t 0 parameterized with t Q . 
They are all on the exact curve W(t) at t = t Q up to 0(e 3 ), but only valid locally for t near 
to- So it is conceivable that the envelope E of {Ct 0 }t 0 which contacts with each local solution 
at t = to will give a global solution. Thus the envelope function x E (t) coincides with W(t); 
x E (t) = x(t,t) = W(t). 

Our task is actually to determine A(to) and 9(to) as functions of t Q so that the family of the local 
solutions has an envelope. According to the classical theory of envelopes given in the previous 
section, the above program can be achieved by imposing that the envelope equation 


dx(t,t 0 ) _ 
dtp 

gives the solution tp = t. From Eq. ’s (E3D, we have 


(3.4) 


dA A d9 e 2 

— + eA = 0, -—I- 

dt 0 dt 0 8 


0, 


(3.5) 


where we have used the fact that dA/dt is 0(e) and neglected the terms of 0(e 3 ). Solving the 
equations, we have A(t 0 ) = Ae~ et °0 and 9(t 0 ) = —+ 9, where A and 9 are constant numbers. 
Thus we get 


x E (t) = x(t,t) = W 0 (t) = A exp(-|t)sin((l - f —)t + 9), 


(3.6) 


up to 0(e 3 ). Noting that y^l — e 2 /4 = 1 — e 2 /8 + 0(e 4 ), one Ends that the resultant envelope 
function x E (t) = W 0 (t) is an approximate but global solution to Ea. dd.llh 
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4 The RG-reduction of dynamics of a generic evolution 
equation 

The theory of the reduction of evolution equations must give a definite method to find out the 
vector field G(s) as well as to construct the attractive manifold M or the function R(s) in Eq. (USD- 
It is often that these tasks can be achieved in a perturbative way as follows, 

d s 

- = G 0 (a) + 1 (a), R(a) = R 0 (s) + p(a), (4.1) 

where Ro(s) gives the coordinate of the unperturbed invariant manifold Mo and p gives the 
deformation of the manifold by the perturbation, as shown in Fig.3. The unperturbed vector 
field Go(s) governs the reduced dynamics on Mo and / y(s) the modification of the dynamics. The 
important point lies in the fact that the modification of the manifold and the dynamics are both 
still function of s parameterizing the unperturbed manifold M 0 . 



Figure 3: The geometrical image of the perturbative construction of the attractive manifold M 
and the reduced dynamics of the vector field G on it. M 0 and Go denote the unperturbed ones. 

Such a view on the reduction of dynamics was emphasized by Kura/moto|14|. It is remarkable that 
Bogoliubov gave the notion of invariant manifold in his contribution to the theory of non-linear 
oscillators [T^; he also described the fluid dynamical limit of Boltzmann equation as a construction 
of a (slow) invariant manifold spanned by the hydrodynamical quantities inbeded in the functional 
space composed of the sigle-particle distribution function [TH], 

In this section, we show how the renormalization-group method works to make the reduction of 
the dynamics of a generic system possessing the possible reduction of dynamics in the perturbative 
waybill: It will be clarified that the system reduction is accomplished by explicitly constructing the 
invariant manifold and the slow dynamics on the manifold in the perturbative way. We emphasize 
that the initial values are chosen by using a simple formula for the special solutions to differential 
equations as in the previous section. 

We treat the following rather generic vector equations in this section: 

d t u = Au + eF(u), (4.2) 

where <9 t u = du/dt, A is a linear operator, F a nonlinear function of u and e is a small parameter 
(| e | < 1). We assume that A has multiply degenerated zero eigenvalues and other eigenvalues of A 
have a negative real part. We assume that A has semi-simple 0 eigenvalues in the present paper; 
the case when A has a multidimensional Jordan cell is also treated in HU- 
We are interested in constructing the attractive manifold M at t —> oo and the reduced dynamics 
on it. We try to construct solve the problem in the perturbation theory by expanding u as 

u(t; t 0 ) = u 0 (t; t 0 ) + eui(t; t 0 ) + e 2 u 2 (f; t 0 ) H-, (4.3) 
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with the initial value W(£ 0 ) at an arbitrary time t 0 . The equations in the first few orders read 


(<% - A)u 0 = 0 , (d t - A)ui = F(u 0 ), (d t -A)u 2 = F'(u 0 )ui, (4.4) 


where 


(F / (u 0 )ui) i = E {<9(F'(uo))i/dMj} (4.5) 

j =i 

is a Freche derivative for the n-dimensional vector. 

We suppose that the equation has been solved up to t — to and the solution has the value W(t 0 ) 
at t 0 . Actually, the initial value must be determined by the perturbative solution self-consistently; 
indeed, u(£) = W (t) is the solution to (14.21) in the global domain. Therefore it should be also 
expanded as follows; W(f 0 ) = W 0 (f 0 ) + e\Vi(t 0 ) + e 2 W 2 (fo) + • • • = W 0 (t 0 ) + p(t 0 ), where p(t 0 ) 
is supposed to be an independent function of Wo. They are not yet known at present but will be 
determined so that the perturbative expansion becomes valid. One of the main purposes in this 
section is to show how sensibly the initial values can be determined order by order. 

In the present paper, we confine ourselves to the case where A has semi-simple 0 eigenvalues. 
Let the dimension of kerA be m; AUj = 0, (i — 1,2,, m). We suppose that other eigenvalues 
have negative real parts; AU Q = Aq,U cv , (a = m + 1 , m + 2, • • •, n), where ReA Q < 0. One may 
assume without loss of generality that Uj’s and U Q ’s are linearly independent. 

The adjoint operator A f has the same eigenvalues as A has; AbU* = 0, (i — 1,2,... ,m) and 
Atu„ = A*U Q , (a = m + 1, m + 2, • • • , n). Here we suppose that TJY and U a ’s are linearly 
independent. Without loss of generality, one can choose the eigenvectors so that (U,;,U Q ) = 0 = 
(U a , Uj), with 1 < i < m and m + 1 < a < n. We denote the projection operators by P and Q 
which projects onto the kernel of A and the space orthogonal to KerA, respectively. 

Since we are interested in the asymptotic state as t —» oo, we may assume that the lowest-order 
initial value belongs to KerA: 


W 0 (t 0 ) = E^o )Ui = W 0 [C}. (4.6) 

i=l 

Thus trivially, u 0 (t; t 0 ) = e^ _ * 0 ^W 0 (t 0 ) = YA!Li Ci(^o)Ui. We notice that a natural parameteriza¬ 
tion of the invariant manifold in the lowest order M 0 is given by the set of the integral constants 
C = C m ) being varied. 

The first order equation (S3) with the initial value Wi(fo) which is not yet determined is formally 
solved to be 


Ul (t;t o) = e^^[W 1 (to) + A- 1 QF(W 0 (t 0 ))] 

+(t ~ t Q )PF(Wo(to)) ~ A- 1 QF(W 0 (t 0 )). (4.7) 

The first term has a possibility to give rise to a fast motion, which should be avoided and are 
analogous to divergent terms in quantum field theory; the divergent terms are subtracted away 
by counter terms, which are analogue to the initial values W, here. Indeed it is nice that the 
initial value Wx(f 0 ) not yet determined can be chosen so as to cancel out the would-be fast term 
as follows; Wi(i 0 ) = —A' 1 QF(W 0 (fo)), which satisfies PWi(t 0 ) = 0 and is a function solely 
of C(t 0 ). Thus we have for the first order solution, Ui(f;t 0 ) = (t — to)PF — A - 1 QF, where the 
argument of F is Wo [C]. Now the invariant manifold is modified to Mi given by 

Mx = {u|u = Wo-eA“ 1 QF(W 0 )}. (4.8) 
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If one stops to this order, the approximate solution reads 


u(t; t 0 ) = W 0 + e{(t - t 0 )PF - A~ l QF}. (4.9) 

Then the RG equation du/dt 0 \ to=t = 0 gives W 0 (t) = ePF(W 0 (t)), which is reduced to an 
m-dimensional coupled equation, 


Ci(t) = e(U t , F(Wo[C])), (i = 1, 2 , • • •, m). (4.10) 

One now sees that ePF(W 0 [C]) gives the vector hied Gq(s) with C being identified with s in 
Ea. (l4.1l) . The global solution representing a trajectory on the invariant manifold up to this order 
is given by 


u(t) = u(t; t 0 = t) = J2 - eA- 1 gF(Wo[C]), (4.11) 

i= 1 


with C(t ) being the solution to (14.101) . 

In short, we have derived the invariant manifold as the initial value represented by (EB and 
the reduced dynamics «nnp on it in the RG method in the first order approximation. 

The above procedure can be easily extended to second and higher orders and the modification 
of the vector held 7(s) in Eq. (l4.1 1) is readily obtained, as shown in [TTj . 

5 Fluid dynamical limit of Boltzmann equation 

In this section, we apply the RG method formulated in the previous sections and in CH to obtain 
the fluid dynamical limit of the Boltzmann equation^]. This is an example of reducing a kinetic 
equation to a slower dynamics]]]?]. In this article, we complete the derivation of the Navier-Stokes 
equation with some corrections to the previous treatment jT2J; we shall also work out to give the 
explicit formulae of the transport coefficients. 

5.1 Basics of the Boltzmann Equation 

The Boltzmann equation^) I2D] is a transport equation which describes the time evolution of 
one-particle distribution function dehned in the phase space: 

| jf( r : v, t) + v ■ V/(r, v, t ) = I[f](r, v, t ). (5.1) 

The right-hand side of the above equation is called the collision integral, 

v , t) = jd 3 v i Jd 3 v 2 Jd 3 v 3 u>(v , t>i|v 2 , v 3 ) (/(r, v 2 , t)f(r , v 3 , t) - f(r , v , t)f(r, iq, i)) ,(5.2) 

where u>(v, V\\v 2 , v 3 ) denotes the transition probability which comes from a microscopic two- 
particle interaction. We remark that the transition probability u>(v, v\\v 2 , v 3 ) contains delta 
functions reflecting energy-momentum conservation law, and satisfies the following relations based 
on the indistinguishability of identical particles and the time reversal symmetry in the scattering 
process: u(v, Ui|u 2 , v 3 ) =a;(«i, v\v 3 , v 2 ) = u>(v 2 , v 3 \v, ui) = u(v 3 , u 2 |ui, v). 

To make explicit the correspondence to the general formulation given in the previous section [TTj. 
one may treat the argument v as a discrete variable [HE]. Discriminating the arguments (r, t) and 





v , we use v as a subscript for the distribution function: /(r, v, £) = /u(r, £) = [/(r, t)} v . Then 
the Boltzmann equation now reads 


t) + » • V/»(r, i) = /[/]»(r, (), 


(5.3) 


where 


I[f]v(r, = <*>(«, v il^ 2 , v 3 ) (fv 2 (r, t)f Vs (r, t ) - /y(r, t)f Vl {r, t)). (5.4) 

Vi v 2 v 3 

As promised, we apply the RG method to extract the low-frequency dynamics from a given 
kinetic equation; in other words, we achieve the coarse-graining of temporal scale by the RG 
method. We introduce e in front of spatial derivative of the Boltzmann equation to make the 
application of the perturbation theory possible: 

| \fv( r , t) = I[f]v(r, t)-ev- Vf v (r, t ). (5.5) 

5.2 Procedure 1 : invariant manifold &; approximate solution 

As was done in the previous sections, we first expand the initial values as follows: 

fv(r, t 0 ) = fv\ r i t 0 ) + ef§\r, t 0 ) + e 2 fg\r , t 0 ) + ■ ■ ■. (5.6) 

Then let jt>(r, t; t 0 ) be an approximate solution around t = t 0 , which obeys Eq. (15. 5 j) with the ini¬ 

tial condition at t = t 0 : /-y (r, to ; t 0 ) = /u(r, t 0 ). We try to solve fv(r, t ; t 0 ) by the perturbation 
theory by expanding it as 

fv{r, t; t 0 ) = f$\r, t; to) +f/?(h t; to) + e 2 /fth t; to) H-, (5.7) 

with the respective initial conditions; fy 1 (r. t 0 ; t 0 ) = fy\r, t 0 ) for n — 0, 1, 2, • • •. Substitut¬ 

ing the above expansion in <ra, we obtain the series of the perturbative equations. The lowest 
few equations read 


d_m 

dt !v 

d ~ ( i) 

dt Iv 

d_ ~( 2) 

9t /v 


= /[/] 




/=/ 


>( 0 ) ’ 


V — 

k df k 

E — 

fc df k 


Ilf] 


Ilf] 


v 


v 


f=f 


f 

(o) J k 


(i) 


/=/ 


f 

7.(0) J k 


( 2 ) 


v ■ V/i 0) , 


® • v/i 11 . 


(5.8) 

(5.9) 
(5.10) 


Here a remark is in order: We have actually used the linearized Boltzmann equation m ne¬ 
glecting the second-order term of / in (15.101) . It is known that the neglected term produces the 
so-called Burnett terms which are absent in the usual Navier-Stokes equations |22j- 


5.3 Procedure 2 : order-by-order analysis 

Here we summarize the significance of each order of equation: 
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• Oth order 

The Oth-order equation is in general a nonlinear algebraic equation which determines the 
lowest-order invariant manifold M 0 (= /^). The special solution of this equation includes 
some integration constants. It turns out that these would-be constants become the slow 
variables on M 0 by the RG equation. 

• 1st order 

The lst-order equation is a linear differential equation because we consider the dynamics 
near Mq. First we derive the zero modes from the eigen vectors of its time evolution operator. 
Then we define the appropriate inner product and projection operator to the kernel space 
spanned by the zero modes. Using these definitions, we solve the linear differential equation. 
The condition that the fast modes orthogonal to the zero modes vanish determine the initial 
condition in this order and thereby the lst-order invariant manifold M\{= /* '*); the lst- 
order perturbation gives the deformation of the invariant manifold from M 0 to HR. The 
coordinates to describe the slow modes are still defined on M 0 . 

• 2nd and higher orders 

The 2nd-order equation is a linear differential equation with the same time evolution oper¬ 
ator. We can solve the 2nd-order equation with the same procedure as the 1st one. The 
2nd-order invariant manifold M 2 (= f ) is determined in the same manner as in the 1st or¬ 
der. This procedure is able to be continued up to arbitrary orders. Although the successive 
deformation of the invariant manifold is obtained, the coordinates for the slow variables are 
still on M 0 . 


Now we are interested in the slow motion which may be realized asymptotically as t —> oo. There¬ 
fore we put the stationary condition on dEIHD and the Oth-order equation becomes the following 
nonlinear algebraic equation: 


9 7(0) 

W v 


m 


v 


f=f 


( 0 ) 


= 0. 


(5.11) 


The Oth-order approximate solution / ( f is the fixed point of the collision integral, and thus is a 
local equilibrium distribution function or Maxwellian: 


fv\r, t; t 0 ) = n(r , t 0 ) 


m 


2vr T(r, t 0 ) 


exp 


m\v — u(r, to 
2 T(r, t 0 ) 


= fv q (r-,to). (5.12) 


Here, the local density n(r, to), local temperature T(r, t 0 ), and local flux u(r, t 0 ) are all dependent 
on the initial time t 0 and the space coordinate r but independent of time f. The Oth-order invariant 
manifold / (0) is given by the initial condition: f$\r, t 0 ) = fv\ r i h); to) = fv l { r 'i to)- The Oth- 
order result is summarized as follows: 


f\t) = n, 

f [0) M = r>. 


The lst-order equation now reads 


d_~ f W 

dt T 


A'f ] + F. 


(5.13) 


(5.14) 
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where the time evolution operator A and the inhomogeneous term F are defined by 


A vk ~ 


d 


d fk 


I[f]v 


f=r 


f v = -v v/: 


eq 

V ! 


(5.15) 


respectively. 

As mentioned above, we must clarify the properties of the linear operator A to proceed further. 
For this purpose, we convert A to the following operator, 

L vk = fv rl A vk fk = ~ £ E E v 3) fv\ ( s vk + S Vl k - S v 2 k - 6 v 3 k)’ ( 5 - 16 ) 

Vt v 2 v 3 


which is called the collision operator [TT! IT^ . 

Let us define the inner product between arbitrary vectors, ip and xf) by 


= Ys fv<Pv (5.17) 

v 

We remark that this definition of the inner product is more adequate than that given in m- a 
nice point is that L becomes self-adjoint with this inner product; (ip, Lx/)) = (L<p,xj>). It is 
essential for the following discussions that L has the zero modes and dim[KerL]= 5[20j: 

Lp>l = 0 , (a = 0,1, 2, 3,4) (5.18) 


where the normalized five vectors <p° a are given as follows with Sv = v — u; 



(5.19) 

(5.20) 

(5.21) 


with ( ip Q a , ip Q p) = Sa/ 3 - The other eigenvalues are found to be negative; in fact one can show that 
( ip , L ip) < 0 for all ip, which means that the kinetic dynamics near the Oth-order solution / eq 
has an attractive slow manifold. 

Next we define the following projection operator P onto the kernel of L , 


P ^\ v = E '*!>), 


(5.22) 


a=0 


and introduce Q = 1 — P as the projection operator to the space complement to KerL. 
Multiplying Eq. dm by the inverse matrix of / eq , we have 

- f rn ) = L (/ «i-i p) + (/«!-> F), 

where /® q . = /u* Efe' l s t-° r( ier solution is readily obtained as 


(5.23) 


/“>(() = \ f m (t 0 ) + A- 1 QF]+(t-t a )PF-A~ 1 QF. 


(5.24) 


where P = f eq P f eq 1 and Q = f eq Q f eq 1 . 
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The first order initial value is now determined so that the would-be fast mode disappear, which in 
turn gives the deformation of the invariant manifold and thereby the lst-order invariant manifold. 
Thus we have for the the lst-order solution 


f<- 1 \to)=-A~ 1 QF, 

/ (1) (t) — (t — t 0 )PF — A~ l Q F. 

Then the 2nd-order equation reads 

Tf (2) = A f <2> + (t-to) H + I, 

where 


H v = — v •V 


PF 


iv 


I v = v • V 


A~' QF 


v 


The solution to this equation is found to be 


;( 2 ) 


f f) = e 


— p (t~to)A 


/ (2) (to) + A~ 2 QH + A" 1 QI 


+ j ~ *o) 2 


PH 


+ {t-to) PI-A^QH - A 2 QH A A~ l QI . (5.28) 


Then the 2nd-order results are summarized as follows: 


P\t) = \(t- i„) 2 


PH 


+ {t-t 0 ) PI — A' 1 Q H - A~ 2 QH + A~ l QI 


(5.25) 


(5.26) 


(5.27) 


(5.29) 


As a result of the above order-by-order analysis, the invariant manifold and the approximate 
solution up to the 2nd order are found to be 


fit 0 ) = f*-eA- x QF-<?[A- 2 QH + A- x Qi 
fit) = r* + e((t-t 0 )PF~A- 1 QF) 


(5.30) 


(t-tof 


PH 


+ (t-t 0 ) PI-A^QH - A~ 2 Q H + A' 1 QI .(5.31) 


Notice the appearance of secular terms in (15.311) . 


5.4 Procedure 3 : envelope equation or RG equation 

Eg. (15. 31 1) shows that the local approximate solution moves away from the invariant manifold as 
\t — to | becomes large owing to the secular terms. The appearance of the secular terms invalidates 
the perturbation expansion of the solution around t ~ f 0 - This evolution is described by the mi¬ 
croscopic time described by the kinetic equation Eq. Q. We can obtain the global solution valid 
in a global domain by constructing the envelope of these diverging local solutions parametrized 
by to- The envelope equation or the RG equation reads: 

= 0, (5.32) 




tn =t 


= 0 OT -£-fv(r, t; to) 
dt 0 


































which is reduced to 


0 = fv ~ e 


PF 


v 


— e 


PI - A' 1 QH 


d j 

- vr:S e 

A" 1 Q F 

+ e 2 

\v dt ( 


v L 


A' 2 QH + A' 1 QI 


-l 


v 


(5.33) 


This RG equation is an equation of motion governing the time evolution of the five slow variables 
n(r, t), T(r, t) and u(r, t ) in / eq . If this equation is solved exactly and these variables are 
substituted in (Enup at to — t, we can construct the macroscopic time evolution of the one- 
particle distribution function f(t). In the discussion below we reduce the master equation (J5.33J1 
to a five-dimensional coupled equation. Applying the projection operator P from the left of (101 
we obtain the following: 


0 = 


Pf 


eq 


V 


— e 


PF 


\v 


— e 


PI 


\v 


[ P^- t (eA- l QF + e 2 (A~ 2 Q H 


A~ l QI 


J v 


(5.34) 


We notice that the time-derivative does hit to the linear operator A which depends on the zero-th 
order distribution function / eq , which case is in contrast to that treated in the previous section; 
we remark that this point was not fully recognized in H3 
Now multipling Eq. (15.341) by the zero modes p° aV an d summing up it in terms of v, we have 

0='EAvfv-c'£,‘e%Fv~e 2 Y. &v h for a = 0, 1, 2, 3, 4. (5.35) 

V V V 

Here we have used the following relations obtained from the definitions (EH and EH- 


E Vav \ P ^] V = E PaV^V, E <plv \Q „ = °> 


V 


and the equality 


\ ^ o ® 

E VaV o 

V Ul 


A~ l QF 


v 


v 


V 


V 


+ e z 


A~ 2 QH + A- l QI\\= 0, 


(5.36) 

(5.37) 


which follows from the fact that <p° a £ Ker L. Notice that (15.351) has the same information as the 
RG equation fl5.33jl . 


5.5 Explicit reduction to Navier-Stokes equation; transport coeffi¬ 
cients 


We now show that Eq. (15.35(1 is nothing but the Navier-Stokes equation |2Dj. Performing the sum¬ 
mation in terms of v, we find that the first and second terms of Ea. (15.351) are evaluated to be 


= (A=—)mn, 

v \y/n m 

fv = (~k= ) mnii 1 for i = 1, 2, 3, 

\y/n V T m' 


(5.38) 

(5.39) 

(5.40) 


and 


E Eo v Pv 

v 

E Fv 

V 

E P4U Fv 

V 


(—= —) m V • (n u), 
\y/n m 


— (—= i/E —) (mnu- Vu 1 + V l (nT)') for i = 1, 2, 3, 
Vyn V T m' V / 


(5.41) 

(5.42) 

(5.43) 
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(5.44) 


In a similar way, the third terms of Eq. (15.3511 are evaluated to be 
E <Pov = 0, 

v 

1 frn 1 

- s/n V T m 


Y.rfvlv = (- r J^-)W{T i *,L- 1 Qf"r 1 F) for * = 1,2,3, 

“ Va/ n V / ra/ 


(5.45) 


V^1V3T 
where T iJ and J' are defined as 


(4= (v l ( r , 27 1 Q r ^- 1 F) + {T\ 77 1 Q r - 1 F ) \ (VV + W)) ,( 5 . 46 ) 


7# = mAffi (Rd, 4 = (4 m |<5u| 2 - ^ t) <5tA (5.47) 

Here we have replaced tp\ v 5v l by J l v because 5v l belongs to KerT and its inner product with 
L~ x Q f “I ’ 1 F 


J v 


is null. 


To proceed further, the explicit representation of L 1 Q f eq 1 F is necessary. By using the 


projection operator lET22t . we have 


fc 


“ C _ “ Sij V • u) + 4 V* In T 


Following the above equation and the space rotational symmetry1171 'ETj[. we arrive at 

( T ij , 77 1 Q / eq “' F) = 2 T) (VV + W - \ V-u), 

(X , 77 1 Q / eq_1 F) = X VT, 

where 77 and A are the so-called transport coefficients defined by 

V = -^T,fvTh 2 L^ k T^ 2 = ^ | o °°dT(T 12 (0),T 12 (r)), [^(r)]^ = E [ f 


vk 


. vk k 


(5.48) 

(5.49) 

(5.50) 

74 (5.51) 


1 1 /*oo r _ r 

A = "T2 E fv J b L vk J k = T 2 L dT ( Jl (°)’ Jl ( r )), [ Ji ( T )]u = ^ [ eiT 

-7 j -*• «/ 0 j 


4- ( 5 - 52 ) 


. ufc ^ 


We note that the transport coefficients, shear viscosity 77 and heat conductivity A, are obtained 
as the correlation function of the microscopic currents (15.4711 like Kubo formula j!8j. 

Putting back e — 1, the reduced RG equations ()5.35|) are found to be 


where 


h + 'V ■ (nu) = 0, 

mnu 1 + mnu ■'Vu 1 = — V J for i = 1 , 2 , 3, 

Tie + 77 it • Ve = - V-J-P ij D ij , 


e = ^ T, p = nT, D ij = ^ ('V V + V J V), 

P ij = S ij p-2r 1 D ij - ( - - 77 ) S ij V • tt, J = -A VT. 

3 


(5.53) 

(5.54) 

(5.55) 


(5.56) 

(5.57) 


These equations (15.5311 - (15.551) is identified with the fluid dynamic equations with dissipation, i.e., 
the Navier-Stokes equation (201 . 

In summary, we have shown that the Navier-Stokes equations can be neatly reproduced as the 
fluid dynamical limit of the Boltzmann equation in the RG method. 
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6 Summary and concluding remarks 

We have shown that the RG gives a powerful and systematic method for the reduction of the 
dynamics and also provides a transparent way for the construction of the attractive slow manifold. 
We have indicated the relation of the underlying mathematics of the RG method with the classical 
theory of envelopes in mathematical analysis. Although the uses of envelopes for physics problem 
was first noted by M. Suzuki in the theory of CAM(Coherent Anomaly Method) [21] for the critical 
phenomena in statistical physics, the usefulness of the notion of envelopes was not well recognized. 

We have obtained the Navier-Stokes equation from the Boltzmann equation by applying the RG 
method: We have worked out for constructing the projection operators and thus explicitly given 
the forms of the transport coefficients in terms of the one-particle distribution function. 

The RG method presented here has a wide range of applicabilities even being confined to the 
transport equations |12| : It can be applied to obtain the Boltzmann equation from the Liouvillc 
equation. The Focker-Planck equation is equally obtained from the Langevin equation by this 
method. Furthermore, the further reduction of the Focker-Planck equation can be also done by 
the RG method. 

As for the reduction of the hydrodynamic equation from the Boltzmann equation, it would 
be interesting to apply the present method to the relativistic case. We also mention that the 
present method is also applicable to extract the critical slow dynamics around the critical point of 
phase transitions [HI EJ. It would be also interesting to apply the method for extracting the slow 
dynamics, say, around the QCD critical end point fZ% 1, 
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